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In order to promote our basic understanding on the Kondo behavior recently observed in europium compounds, we 
analyze an impurity Anderson model with seven / electrons at an impurity site by employing a numerical renormaliza¬ 
tion group method. The local part of the model consists of Coulomb interactions among / electrons, spin-orbit coupling 
A, and crystalline electric field (CEF) potentials, while we consider the hybridization V between local / electrons and 
single-band conduction electrons with symmetry. For A = 0, we observe the underscreening Kondo behavior for 
appropriate values of V, characterized by the entropy change from In 8 to In 7, in which one of seven / electrons is 
screened by conduction electrons. When A is increased, we obtain two types of behavior depending on the values of 
V. For large V, we find the entropy release of In 7 at low temperatures, determined by the level splitting energy due to 
the hybridization. For small V, we also observe the entropy change from In 8 to In 2 by the level splitting due to the 
hybridization, but at low temperatures. In 2 entropy is found to be released, leading to the Kondo effect. We emphasize 
that the Kondo behavior for small V is observed for realistic values of A in the order of 0.1 eV. We also discuss the effect 
of CEF potentials and the multipole properties in the Kondo behavior found in this paper. 


1. Introduction 

Research on heavy-electron systems has been one of central 
topics in the field of strongly correlated electron physics. 
The origin of such heavy-electron state has been understood 
on the basis of quantum critical phenomena,emerging 
in the competing region of itinerant properties of / electrons 
due to the Kondo effect"^' and localized nature due to the 
Ruderman-KitteTKasuya-Yosida (RKKY) interaction. 

The Kondo effect has been first discovered as resistance 
minimum phenomena in good metals such as Cu, Ag, and 
Au with a small amount of magnetic impurities such as Mn 
and Fe. In the dilute system, it has been widely recognized 
that Kondo effect occurs when the singlet state is formed 
from local magnetic moment due to the antiferromagnetic 
(AF) coupling with conduction electrons.As for the orig¬ 
inal problem of resistance minimum phenomena in metals 
with magnetic impurities, Kondo has actually demonstrated it 
by quantum-mechanical calculations for scattering amplitude 
of electrons due to magnetic impurities. 

On the other hand, there appears long-range interaction be¬ 
tween localized electrons via conduction electrons. This is 
called the RKKY interaction,^ which enhances the local¬ 
ized nature of / electrons. As a result of the competition be¬ 
tween the Kondo effect and the RKKY interactions, quantum 
criticality appears between AF and metallic phases, when we 
control the coupling constant pJcf, where Jet denotes the AF 
coupling between localized and conduction electrons and p 
is the density of states at the Fermi level. The picture has 
been summarized in the Doniach’s phase diagram,'^' which 
has been a guiding principle in the heavy-electron physics to 
discover the unconventional superconductivity mediated by 
quantum critical fluctuations. 

Along this line, the heavy-electron phenomena and uncon¬ 
ventional superconductivity have been intensively and exten¬ 
sively investigated in the Ce compounds with one / electron 
for Ce^+ ion, since the pioneering discovery of supercon¬ 
ductivity in CeCu 2 Si 2 .*®^ Recently, the research on heavy- 


been activated on the basis of the electron-hole conversion 
picture of p has been claimed that this material 

exists just on the quantum critical point at ambient pressure. 
For the case of two / electrons, there is a long history in 
the research of U compounds. In fact, after the discovery 
of heavy-electron superconductivity in CeCu 2 Si 2 , supercon¬ 
ductivity has been found in U compounds such as UBeia,^^' 
UPt3,24> URu2Si2,2^^ UPd2Al3,2®) and UNi2Al3.2^> In a re¬ 
cent decade, superconductivity and magnetism in Pr com¬ 
pounds such as PrOs2Sbi2,2*^) PrPb3,2‘^’^°> and PrT 2 X 2 o (T: 
transition metal, X=A1 and have been actively in¬ 

vestigated. On the basis of the electron-hole conversion pic¬ 
ture of Pr^+, it is also interesting to note Tm^+ with twelve / 
electrons. In TmsRhgSnig, superconductivity has been found 
with Tc = 2.2K.^^' Peculiar reentrant properties have been 
considered to be related to the coexistence of magnetism and 
superconductivity. 

When we look back over the history of heavy-electron ma¬ 
terials, we immediately notice that elements around both ends 
of lanthanide series have been focused thus far. However, 
recently, the heavy-electron state in Eu compounds, just at 
the center in lanthanide series, have attracted renewed atten¬ 
tion.In the divalent ion of europium, seven electrons are 
accommodated in the 4/ orbitals, while in the trivalent ion 
of europium, we find six 4/ electrons. In an LS coupling 
scheme, due to the Hund’s rules, we obtain J = S = 7/2 
and L = 0 for Eu^+, where J, S, and L indicate, respec¬ 
tively, total angular momentum, total spin momentum, and 
total orbital momentum. On the other hand, we find J = 0 
with S = L = 3 for Eu^+, which is non-magnetic. The dif¬ 
ference in magnetic properties between Eu^+ and Eu^+ ions 
are significant, but in general, the valence fluctuations easily 
occur in Eu compounds, since the energy difference between 
two valence states has been known to be small. Thus, the va¬ 
lence fluctuation is one of key issues in Eu compounds. 

Quite recently, in Eu 2 Ni 3 Ge 5 and EuRhSi 3 , heavy-electron 
states have been claimed to be observed in the temperature 
denendence in the measurement of resistivitvf^^l At ambient 
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(b) LS coupling scheme 



Fig. 1. Electron configurations of the seven /-electron state for (a) the j-j 
coupling scheme and (h) the LS coupling scheme. Note that up and down 
arrows denote, respectively, pseudo-spin up and down electrons in (a), while 
in (h), they indicate real-spin up and down electrons. The pseudo-spin state 
is defined through the time-reversal relation. 


pressure, those compounds are in the AF state at low temper¬ 
atures, but under pressure such as several GPa, the AF state 
is suppressed and three characteristic temperatures have been 
observed in the resistivity. The highest one was assigned as 
Tv, which is the valence transition temperature. Other two low 
temperatures were considered as Kondo temperatures. These 
results were claimed to be quite similar to those of Ce com¬ 
pounds. 

At a first glance, it seems to be difficult to accept the simi¬ 
larity between Ce and Eu compounds, but it is easy to hit upon 
an idea on the basis of a j-j coupling scheme, where j denotes 
the total angular momentum of one / electron. Namely, in the 
j-j coupling scheme, among seven / electrons, six electrons 
fully occupy the sextet of j = 5/2, whereas one electron is 
accommodated in the octet of j = 7/2, as shown in Fig. 1(a). 
Since six electrons in the sextet do not contribute to electronic 
properties, one / electron plays a main role for the Kondo ef¬ 
fect, leading to the similar behavior as that of Ce compound. 
This idea has been also emphasized by the present author for 
the explanation of active quadrupole degrees of freedom in 
Gd compound with seven / electrons in the trivalent ion state 
of gadolinium. 

Here we cast a naive question: Are there any problems to 
use the j-j coupling scheme for the Eu compounds? In fact, 
in a textbook of solid state physics, it is standard to use the 
LS coupling scheme. Namely, first we construct the many / 
electron state characterized by S and L due to the so-called 
Hund’s rules. Then, we consider the effect of spin-orbit cou¬ 
pling by including the term of KL ■ S with a spin-orbit cou¬ 
pling A, leading to the multiplet characterized by J. Eor seven 
/ electrons, as mentioned above, first we obtain the state with 
5 = 7/2 and L = 0 from the Hund’s rules, as shown in 
Eig. 1(b). Since L is zero, the ground state is characterized by 
J = S = 7/2. 

Eor rare-earth compounds, in general, the magnitude of 
Hund’s rule interaction among / orbitals is a few eV, while the 
spin-orbit coupling takes a value of a few thousand Kelvins. 
If we should take one of the two limiting situations, it is better 
to choose the LS coupling scheme, as readers have learned in 
the standard textbook. In this sense, it seems to be difficult to 
understand the Kondo effect in the Eu compound on the basis 
of the j-j coupling scheme. 

However, we should remark that both of the spin-orbit cou¬ 
pling and the Hund’s rule interaction are finite in actual ma¬ 
terials and the actual situation is always in the middle of the 


LS and j-j coupling schemes. Namely, the wave function of 
the many / electron state is the mixture of those in the LS 
and j-j coupling schemes. In order to discuss the Kondo ef¬ 
fect in the Eu compound even qualitatively, it is essential to 
consider both of the spin-orbit coupling and the Hund’s rule 
interaction. This point has been emphasized in the research of 
quadrupole susceptibility in Gd compounds. 

In this paper, we analyze the seven-orbital Anderson model 
by employing a numerical renormalization group technique. 
The local term contains Coulomb interaction, spin-orbit cou¬ 
pling, and crystalline electric field (CEE) potential terms. 
Here we introduce the Hund’s rule interaction Lf, the spin- 
orbit coupling A, and the CEE potential W. As for conduc¬ 
tion electrons, we include only one conduction band with au 
symmetry. Readers may consider that a trivial result of the un¬ 
derscreening Kondo effect is obtained, but such a well-known 
result is observed only in the LS coupling limit. When we 
increase the spin-orbit coupling in this situation, we confirm 
the Kondo bahavior similar to that in the Ce compound, as 
expected from the j-j coupling scheme. An important find¬ 
ing in this paper is that the Kondo effect similar to the case 
of n = 1 occurs for a realistic value of spin-orbit coupling 
even in the case of n = 7, where n denotes local /-electron 
number. When we explicitly include the cubic CEE potentials, 
we discuss the multipole properties and the Kondo behavior 
similar to Ce compounds. 

The organization of this paper is as follows. In Sec. 2, 
we provide the local /-electron Hamiltonian and briefly re¬ 
view the change of the seven /-electron states between the 
LS and j-j coupling schemes by evaluating the Curie con¬ 
stant. We emphasize that the transition region between them 
just corresponds to the situation of actual /-electron mate¬ 
rials. In Sec. 3, we introduce the impurity Anderson model 
to discuss the Kondo phenomenon for the case of n = 7. 
We also briefly explain the method used in this paper and 
provide the definition of multipole susceptibilities. In Sec. 4, 
we exhibit our numerical results and discuss how an entropy 
changes with the decrease of a temperature. Eirst we clearly 
show the underscreening Kondo effect in the LS coupling 
scheme for X = W = 0. Then, we investigate the effect of 
the spin-orbit coupling on the underscreening Kondo behavior 
for W = 0. Eor small V and A/C/ in the order of 0.1, we con¬ 
firm the Kondo behavior characterized by the entropy release 
of In 2, as found in the Ce compound. Then, we also investi¬ 
gate the Kondo behavior for W ^ 0 and discuss the multipole 
susceptibility to confirm that the relevant multipole is dipole 
in the present Kondo effect. Einally, in Sec. 5, we provide a 
few comments on future problems and summarize this paper. 
Throughout this paper, we use such units as = fee = 1- 


2. Local f Electron State 

2.1 Local Hamiltonian 

Eirst, we define the local /-electron Hamiltonian as 
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where fma- denotes the annihilation operator for local / elec¬ 
tron with the spin a and z-component m of the angular 
momentum £ = 3, a = -fl (—1) for up (down) spin, 
fmim 2 ,m 3 m 4 indicates the Coulomb interaction, A is the spin- 
orbit coupling, and Bm,m' denotes the CEF potential. 

The Coulomb interaction I is known to be expressed as 

6 

.^mi m 2 ,77137724 — E F'^Ck{mi,m4)ck{m2,m3), ( 2 ) 

k=0 

where indicates the Slater-Condon parameter and is 
the Gaunt coefficient."^^* Note that the sum is limited by the 
Wigner-Eckart theorem to fc = 0, 2, 4, and 6 . Although the 
Slater-Condon parameters should be determined for the mate¬ 
rial from the experimental results, we assume the ratio among 
the Slater-Condon parameters as 

= 1017, = 517, F^ = 317, F® = 17, (3) 

where U is the Hund’s rule interaction among / orbitals. 

Each element of for the spin-orbit coupling is given by 

Cm,fT;m,(T 777.(7/2, 

,_ (4) 

Cm-t-cr, —(T;m,(T V 1 ) 777 ( 777 -f Cf') j 2,^ 

and zero for other cases. The CEF potentials for / electrons 
from the ligand ions is given in the table of Hutchings for 
the angular momentum £ = 3."*^* For cubic structure with Oh 
symmetry, Bm,m' is expressed using a couple of CEF param¬ 
eters, F® and Fg, as 

Fa,3 = F_3,-3 = 180F2 + ISOF®, 

F 2,2 = F_ 2,-2 = -420F2 - 1080F®, 

Fi 1 = F_i _i = 60F2 -b 2700Fg, 

(5) 

Fo.o = 360F^ - 3600Fg, 

Fa.-i = F_3 .i = 60yT5(F2 - 21F0), 

F2,_2 = 300F2 + 7560FO, 

Note the relation Bm,m> = Bm',m- Following the traditional 
notation,"*®* we define F® and Fg as 

F® = Wx/Fi4), F° = W(1 - \x\)/F{6), ( 6 ) 

where x and y specify the CEF scheme for the Oh point group, 
while W determines the energy scale for the CEF potential. 
We choose F(4) = 15 and F( 6 ) = 180 for £ = 3."*®* 

Before proceeding to the results, we provide a comment on 
the energy scale for U, A, and W, Among them, the largest 
one is U and its magnitude is in the order of 1 eV, since it de¬ 
notes the Hund’s rule interaction among / orbitals. The next 
one is A and its magnitude is in the order of 0.1 eV, although 
the precise values depend on the kind of lanthanide and ac¬ 
tinide ions between 0.1 and 0.3 eV. The smallest one is W, 
since the magnitude is typically in the order of meV, although 
the values depend on materials. In any case, it is reasonable to 
consider F > A ^ IE and A/F ^ 0.1 for actual /-electron 
compounds. 

Concerning the energy unit, we define it as a half of the 
conduction bandwidth in this paper, as will be mentioned 
later. Throughout this paper, we set it as 1 eV, but this value 
is realistic as a half of the conduction bandwidth. 



Fig. 2. Curie constant C for n = 7 as a function of A/f/ for the case of 
W = 0. We find the transition from the value for the LS coupling limit 
to that for the j-j coupling limit in the region pfO.l < X/U < 1, which 
includes the parameters for actual materials. 


2.2 LS and j-j coupling schemes 

In the standard textbook for solid state physics, it is rec¬ 
ommended to employ an LS coupling scheme for /-electron 
states in the limit of F ^ A, while we should use a j -j cou¬ 
pling scheme in the limit of F <C A. Since the actual situation 
is found for A/F ^ 0.1, it seems to be better to choose al¬ 
ways the LS coupling scheme. However, as has been empha¬ 
sized in our previous papers,"*^’the wavefunction of the 
/-electron state for A/F ~ 0.1 is well approximated by that 
in the j-j coupling scheme. The /-electron state is continu¬ 
ously changed from the LS to the j-j coupling limits, when 
we change the ratio of A/F, but the transition region is found 
in the range between A/F = 0.1 and 1. 

For /"^-electron case, we also find such a transition in the 
same parameter region of A/F. In order to reconfirm this 
point for the case of n = 7, we evaluate the Curie constant 
C on the basis of Fioc for W = 0. We obtain the magnetic 
susceptibility x as x = C/T and the Curie constant C is given 
by C = (pj/j,B)^J(7-l-l)/3fcB- Here pb indicates the Bohr’s 
magneton, J denotes the size of total angular momentum J, 
which is given by J = L -\- S with total angular momentum 
L and total spin momentum S and gj is the Fande’s (/-factor. 
Note that this expression for C is in common with the LS and 
j-j coupling limits, although the magnitude of C is changed. 

In the LS coupling limit, for the case of 77 = 7, we obtain 
J = S = 712 and L = 0 due to the Hund’s rules. In this case, 
we find gj = 2, which is just equal to the electron’s g factor. 
Thus, we obtain the magnetic moment as M = 7pB and C = 
21/7|/fcB- On the other hand, in the j-j coupling limit, we 
accommodate one / electron in the octet j = 7/2, while the 
sextet j = 5/2 is fully occupied. In this case, the total angular 
momentum J is, of course, 7/2, which is originating from one 
/ electron in j = 7/2 octet. The Fande’ p-factor of j = 7/2 
is equal to gj = 8/7 and thus, we obtain M = AjiB and C = 
(48/7)/7|/fcB- Due to the reduction of the magnetization, C 
is changed from 21 to 48/7. 

Concerning the value of C between the LS and j-j cou¬ 
pling limits, it is necessary to perform the numerical evalua¬ 
tion of eigenenergies and eigenstates of Fioc for 77 = 7 and 
IE = 0 by changing the ratio of A/F. The result is summa- 
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Fig. 3. (a) Eigenenergies vs. x for A = 0, C/ = 1, and W = —10“^. (b) 

Eigenenergies vs. x for A = 0.2, U = 1, and W = —10“^. (c) Eigenener¬ 
gies vs. a: for A = 10“^, U = 1, and W = —10“^. Note that (a) denotes the 
result in the LS coupling scheme, while (c) indicates that in the j-j coupling 
limit. 


rized in Fig. 2. As mentioned above, we find (7 = 21 and 48/7 
in the limit of X/U = 0 and X/U = oo, respectively, lead¬ 
ing to the good check of the numerical calculations. We note 
the continuous decrease in C from 21 to 48/7, which orig¬ 
inates from the reduction of the magnetic moment M from 
M = 7/rB in the LS coupling limit to M = 4fiB in the j-j 
coupling limit. 

It is emphasized here that the reduction of C occurs in the 
range of 0.1 < X/U < 1, which just corresponds to the values 
of actual /-electron materials. In this transition region, the /- 
electron wavefunction is given by the mixture of those in the 
LS and j-j coupling limits. Thus, the properties of both limits 
can be observed in the range of 0.1 < X/U < 1, leading to the 
expectation of a possible explanation of the Kondo behavior 
in electron systems. 

2.3 CEF level schemes 

Another evidence for the mixture of both limits is found in 
the CEF level schemes. In Figs. 3, we summarize the varia¬ 


tion of the CEF level schemes, when A is increased for the 
fixed values of (7 = 1 and W = —10“^. Note that Fig. 3(a) 
indicates the result in the LS coupling limit, in which we do 
not find any effects of CEF potentials on the /-electron low- 
energy states for n = 7. Readers may be doubtful of this re¬ 
sult, but we should note that J is purely given by spin, i.e., 
J = S' = 7/2 in the LS coupling limit for n = 7. When 
we recall the fact that the CEF potential acts on charge, not 
on spin, there occurs no effect of CEF potentials on the states 
with J = S = 7/2, except for the shift of the total energy. 
Thus, we still obtain the octet independent of CEF potentials. 

Next we consider the j-j coupling limit. Here we pay our 
attention to Fig. 3(c) by skipping Fig. 3(b). In Fig. 3(c), we 
exhibit the result for A = 10^. This is, of course, unphysical 
value, but we use it for the purpose to realize the j-j coupling 
limit. As mentioned above, in this limit, we accommodate six 
electrons in the sextet of / = 5/2 and one electron in the octet 
of / = 7/2. In sharp contrast to the LS coupling limit with 
J = S = 7/2, one / electron in the / = 7/2 octet feels 
the CEF potentials, as shown in Fig. 3(c). Under the cubic 
CEF potentials, the octet is split into two doublets (F/" and 
F/) and one quartet (F^). It is quite natural that the present 
level scheme is quite similar to that for J = 7/2‘^^'* due to the 
symmetry reason, although there are small deviations due to 
the difference in the definition of the CEF parameters. 

Let us turn our attention back to Eig. 3(b). This is the result 
for A = 0.2, which is near the realistic values for Eu and 
Gd ions. Eirst we notice that the energy scale of Eig. 3(b) is 
smaller in one order than that of Eig. 3(c). The CEE energy 
splitting is in the order of meV, which is about 10 K. This 
value is considered to be small, but the CEE excitation in this 
order can be detected in the experiment. Thus, as mentioned 
in the previous paper,it is a challenging issue to measure 
the CEE excitation in Gd and Eu compounds, which have not 
been expected to show CEE excitations. 

Second the symmetries of the CEE ground states are the 
same as those in the j-j coupling scheme. It is quite natural 
that the ground-state multiplet is always characterized by the 
total angular momentum J = 7/2, irrespective of the values 
of U and A. Note that the magnitude of the matrix elements 
among the states of is affected by U and A, where Jz de¬ 
notes the z-component of J. 

3. Model and Method 

3.1 Seven-orbital Anderson model 

Now we consider the conduction electron hybridized with 
localized electrons. In general, all / orbitals are hybridized 
with conduction electrons with the same symmetry, but it is 
very difficult to perform numerical calculations by taking into 
account seven conduction bands. Since the purpose here is 
to reveal the effect of spin-orbit coupling on the Kondo phe¬ 
nomena in /^ electron systems, it is reasonable to consider 
the minimum model to investigate the Kondo phenomena. In 
this sense, the minimum model should include one conduc¬ 
tion band, since we can expect the appearance of the under¬ 
screening Kondo phenomena even for one conduction band 
hybridized with local / orbitals. Then, it is natural to consider 
Ou conduction band, since local Ou state is non-degenerate 
even under the high symmetry ligand field such as the cu¬ 
bic CEE potential. Note that the local Ou state is described as 

/ pf . ... 
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Then, the model is expressed as 

H = E h.c.) + Tfioc, (7) 

k,(T k,m,(7 

where Sk is the dispersion of conduction electron, Cka is an 
annihilation operator of conduction electron with momentum 
k and spin a, and Vm is the hybridization between conduction 
and localized electrons. Note that Vm is given by 14 n =2 = 
— Vm=-2 = V and zero for other components. The energy 
unit is a half of the conduction bandwidth, which is set as 1 
eV throughout this paper, as mentioned above. 

3.2 Numerical renormalization group method 

For the diagonalization of the impurity Anderson model, 
we employ a numerical renormalization group (NRG) 
method,^^’"^*^ in which we logarithmically discretize the mo¬ 
mentum space so as to include efficiently the conduction elec¬ 
trons near the Fermi energy. The conduction electron states 
are characterized by “shell” labeled by N and the shell of 
TV = 0 denotes an impurity site described by the local Hamil¬ 
tonian. Then, after some algebraic calculations, the Hamilto¬ 
nian is transformed into the recursion form as 

Hn+1 = VAHn + Ctv ^(c]v^CAr+lcr + (8) 

cr 

where A is a parameter for logarithmic discretization, cnct 
denotes the annihilation operator of conduction electron in the 
A^-shell, and indicates “hopping” of electron between N- 
and (TV + l)-shells, expressed by 

^ ^ (1+A-1)(1-A-a^-1) 

2y'(l-A-2^-1)(1- A-2W-3)' 

The initial term Hq is given by 

Ho = ^ Vmicljma + fLcO.)]- (10) 

mcT 

For the calculations of thermodynamic quantities. We eval¬ 
uate the free energy F for local / electron in each step by 

F = -r(lnTre-^'^/'^ - InTre”^^/'^), (11) 

where a temperature T is defined as T = in the 

NRG calculation and denotes the Hamiltonian without 
the hybridization term and iVioc- Then, we obtain the entropy 
5'imp by 'Simp = —OF/dT and the specific heat Cimp is eval¬ 
uated by Cimp = —Td'^F/dT'^. In the NRG calculation, we 
keep M low-energy states for each renormalization step. In 
this paper, we set A = 5 and we keep M = A, 500 low-energy 
states for each renormalization step. 

3.3 Multipole susceptibility 

For the purpose to discuss the multipole properties later, 
we provide the definition of multipole operator in this subsec¬ 
tion. Since we discuss the effect of the spin-orbit coupling for 
\/U = 0 to A/C/ = oo, it is necessary to use the same defi¬ 
nition both in the LS and j-j coupling schemes, irrespective 
of the values of A/C/. Thus, we define the multipole as a spin- 
orbital density in the form of a one-body operator from the 
viewpoint of the multipole expansion of electron density in 
electromagnetism. On the basis of this definition of the mul¬ 
tipole operator, we have developed microscopic theories for 

_____ _.t__ __.t__i_ __i_ r 


in Np 02 has been clarified by evaluating multipole interac¬ 
tion by the standard perturbation method in terms of electron 
hopping.We have discussed possible multipole states of 
filled skutterudites by analyzing the multipole susceptibility 
of a multiorbital Anderson model based on the j-j coupling 
scheme.We have also discussed the multipole state in 
actinide dioxides^®'and Yb compounds.®'^ Recently, a mi¬ 
croscopic theory for multipole ordering from an itinerant pic¬ 
ture has been developed on the basis of a seven-orbital Hub¬ 
bard model with spin-orbit coupling.®^' 

The multipole operator T is expressed in the second- 
quantized form as 


FW _ \ ' r^{k) rji{k.q) A £ 

^i,7 / . ^-f,q-^7na,7n'a' J imaMru'cr' ^ 

q, 171 ( 7 ,m'<T' 


( 12 ) 


where k indicates the rank of the multipole, q denotes an in¬ 
teger between —k and fc, 7 is a label used to express an Oh 

(k) 

irreducible representation, is the transformation matrix 
between spherical and cubic harmonics, and can be 

calculated, using the Wigner-Eckart theorem, as®^' 


rp{k,q) 

171 ( 7 ,m' (T' 


E 


•\/2/ + 1 


{j ii\j p.'kq) 




(13) 


Here, f = 3, s = 1/2, j = £ ± s, /r denotes the z-component 
of j, {jp\j'p' j"p") indicates the Clebsch-Gordan coefficient, 
and (/I 11 /) is the reduced matrix element for a spherical 

tensor operator and is given by 




y / (2j + fc + i)! 

2’^^ {2j-k)\ 


(14) 


Note that k < 2j and the highest rank is 2j. Thus, we treat 
multipoles up to rank 7 for / electrons in this definition. 

Here we should note that multipoles belonging to the same 
symmetry are mixed in general, even if the rank is different. 
Namely, the /-electron spin-charge density should be given 
by the appropriate superposition of multipoles, expressed as 


(15) 

fe,7 


where we redefine each multipole operator so as to satisfy the 
orthonormal condition of^b 


Tr{fWf^f')} (16) 

In order to determine the coefficient p^\ it is necessary to 
evaluate the multipole susceptibility in the linear response 

(k) 

theory. Namely, pj is determined by the eigenstate of the 
largest eigenvalue of susceptibility matrix, given by 


Xk-i,k'-f' —■ 


1 ^ - e-^AT 

Z ^ 

1,3 


Ej - Ei 


(*|[f«-pW]|j) 


(17) 




where Z is the partition function given by Z='^^ e ®«/^, Ei 
denotes the eigenenergy of the *-th eigenstate |j) of the Hamil¬ 
tonian, T is a temperature, and {i\Tj^^ |i) /Z. 


nru,^ -^,,14.:_- -i _...:-1— ,.4:4.1.,^ 
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susceptibility matrix. By using the NRG method, we evalu¬ 
ate the matrix elements of the multipole susceptibility in each 
NRG step. Then, we can determined the optimized multipole 
state in an unbiased manner. 

Note that in this paper, we express the irreducible represen¬ 
tation of the CEF state by Bethe notation, whereas for multi¬ 
poles, we use short-hand notations by the combination of the 
number of irreducible representation and the parity of time re¬ 
versal symmetry, “g” for gerade and “u” for ungerade. For Oh 
symmetry, we have nine independent multipole components 
as Ig, 2g, 2u, 3g, 3u, 4g, 4u, 5g, and 5u. Note that lu does not 
appear within the rank 7. 

4. Calculation Results 

4.1 Underscreening Kondo effect 

Fet us show our numerical results. First we consider the 
situation ofU^O and X = W = 0, which is the limit of LS 
coupling scheme. In this case, the ground-state multiplet is 
characterized by S' = 7/2 and L = 0. When one conduction 
band is hybridized with such localized state, it has been well 
known that the underscreening Kondo effect occurs. Namely, 
one of seven electron among the octet forms the singlet due to 
the hybridization with single conduction electron band, while 
another six electrons are still localized. Thus, even after the 
underscreening Kondo effect, the septet remains. Here we use 
the Kondo temperature, even when the underscreening Kondo 
effect occurs. 

In Fig. 4(a), we show the temperature dependence of en¬ 
tropy for several values of V. For F=0.5, we find the entropy 
In 8 in the present temperature range. If we decrease the tem¬ 
perature, we expect the change of the entropy, but it is diffi¬ 
cult to obtain the results with a reliable precision. Then, we 
increase the value of V to elevate the Kondo temperature. As 
expected, we find the gradual decrease in the entropy when 

V is increased. For V=2.0, the entropy immediately becomes 
In 7 in the present temperature range, since the Kondo tem¬ 
perature is high enough. 

Although it is difficult to define the Kondo temperature 
only from Fig. 4(a), it is possible to define Tk as a peak 
temperature in the specific heat, which is formed by the en¬ 
tropy change from In 8 to In 7. The results are shown in 
Fig. 4(b). Note that the absolute values of the specific heat 
are suppressed, since the entropy release is relatively small in 
this case. Tiny fluctuations in the specific heat are observed, 
since the magnitude of the specific heat is relatively small in 
this case and the effect of numerical precision appears. For 

V = 0.5, we find Cimp is zero and no peak is found in the re¬ 
gion of T < 10“^. When V is increased, the value of Cimp is 
gradually enhanced. For V = 0.8 and 0.9, we find the broad 
peak structure in the specific heat. These peaks are consid¬ 
ered to originate from the entropy change from In 8 to In 7. 
For V = 1.0, we observe a weak shoulder structure around at 
T = 10“^, but the peak structure is smeared. For V = 2.0, 
Cimp becomes immediately zero and no peak structure can be 
found for T < 10“^. 

In Oder to confirm that the entropy change and the peak 
formation in the specific heat are indications of the Kondo ef¬ 
fect, we plot the peak temperature Tp in Fig. 4(c). We can 
define the peaks in the region of 0.75 < V < 1, but we 
find that Tp is well fitted by bexp{—a/V^) with appropriate 



Fig. 4. (Color online) (a) Entropy Simp temperature T for U = 1 and 
A = VC = 0. The hybridization V is changed between 0.5 and 2.0. (b) 
Specific heat Cimp vs. temperature for the same parameters in (a), (c) The 
peak temperature Tp (solid symbols) vs. V. The red curve denotes InTp = 
b — ajV^ with appropriate constants a and b. 


mula for the Kondo temperature Tk = Dexp[—l/{pJc{)], 
where D is the half of the conduction electron band. Note 
that D is set as unity in this paper. In the single-band An¬ 
derson model, we can obtain Jcf = 4V^//, where I de¬ 
notes the on-site Coulomb interaction. In the present case, it 
is difficult to derive such a simple form of Jcf, since we con¬ 
sider the complicated orbital-dependent interactions through 
the Slater-Condon parameters. However, we can deduce that 
Jcf can be obtained, in any case, in the second-order pertur¬ 
bation in terms of V. We expect that Tp is in proportion to 
exp(—a/V^) with an appropriate constant a. We also intro¬ 
duce another constant b to fit Tp as 6exp(—a/V^). From the 
result in Fig. 4(c), we conclude that the underscreening Kondo 
effect actually occurs and it is characterized by the entropy 
change from In 8 to In 7 due to the screening of one spin by 
single conduction band. 
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Fig. 5. (Color online) (a) Entropy Simp (solid square) and specific heat 
Cimp (solid triangle) for A = 0.1, U = 1, W = 0, and V = 3.0. (b) The 
peak temperature Tpi (solid circle) and the excitation energy Ai ?2 (solid 
line) of the two-site system as functions of A. Inset shows the schematic view 
for the energy levels of the two-site system, in which the septet is split into 
one singlet and two triplets. 


4.2 Effect of spin-orbit coupling 

Next we include the spin-orbit coupling. Namely, we con¬ 
sider the finite values of A for {7 = 1 and W = 0. Depend¬ 
ing on the values of V, we consider two typical situations for 
large and small values of V such as V = 3.0 and 0.6. From 
Fig. 4(a), the underscreening Kondo effect occurs in the high- 
temperature region as T > 0.1 for V > 2.0, while the octet 
still remains in the temperature range for V < 0.6. 

In Fig. 5(a), we show the typical results for entropy and 
specific heat for A = 0.1, {7 = 1, and W = 0 with the large 
value of V such as = 3.0. As expected from the results in 
Figs. 4, we observe the plateau of In 7 for T < 0.1. When 
we further decrease the temperature, we find the release of 
entropy In 7, leading to a large peak in the specific heat, as 
shown in Fig. 5(a). In the present case, we find such a peak at 
T K. 10“^. Here we define the peak temperature as Tpi. We 
are interested in the origin of this entropy change, but it seems 
to be irrelevant to the Kondo effect. Intuitively, the entropy 
change in this case is rather rapid in comparison with the that 
in the Kondo effect. A more concrete discussion can be done 
by the direct comparison of the relevant energy scales. In one 
word, this is the level splitting due to the hybridization. 

At present we do not explicitly include the CEF potential, 
but it has been well known that the level splitting occurs due 
to the effect of the hybridization between localized and con- 
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Fig. 6. (Color online) (a) Entropy Simp (solid square) and specific heat 
Cimp for A = 0.2, U = 1, W = 0, and V = 0.6. (b) The peak temperature 
Tp 2 (solid circle), the Kondo temperature Tk (solid square with solid line), 
and the excitation energy (broken curve) of the three-site system as 

functions of A. Inset shows the schematic view of the energy levels of the 
three-site system, in which the octet is split into one quartet and two doublets. 


perature Tpi is well fitted by AE 2 , where AE 2 denotes the 
excitation energy between the ground state singlet and the 
first excited triplet in a two-site system, as shown in the in¬ 
set of Fig. 5(b). Note that the two-site system is composed of 
impurity site and one conduction site, which are connected by 
the hybridization V. For the case of A = 0, we find that the 
ground state is septet. Due to the combination of S' = 7/2 
at impurity site and S = 1/2 at the conduction site, we ob¬ 
tain Stot = 4 and 3, but the state with Stot = 3 becomes 
the ground state with septet. Here Stot denotes the magni¬ 
tude of total spin moment. When we include the effect of A, 
the level splitting due to the hybridization becomes significant 
and AE 2 is increased with the increase of A. The temperature 
to characterize the entropy release of In 7 is deduced to be 
determined by AE 2 , as shown in Fig. 5(b). 

Next we consider the case with small V. In Fig. 6(a), we de¬ 
pict the entropy and specific heat for A = 0.2, U = 1,W = 0, 
and V = 0.6. As expected from the discussion in Figs. 5, at 
high temperatures as 10“^ < T < 10“^, we find a plateau of 
In 8 originating from the local octet of J = 7/2. For the case 
of A = 0, only one electron spin is screened by one conduc¬ 
tion band, leading to the entropy change of In 8 to In 7, but in 
the present case with A = 0.2, first we find the entropy change 
from In 8 to In 2 by the level splitting due to the hybridization 
with conduction electrons. Then, the remained In 2 from the 
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local doublet is eventually released by the Kondo effect. 

We claim that the result in Fig. 6(a) is important. The un¬ 
derscreening Kondo effect in the /^-electron system has been 
well understood and probably, it has not been considered to be 
an intriguing phenomenon. However, under the effect of the 
spin-orbit coupling, the entropy of In 8 from the local octet 
due toJ = S' = 7/2is released in the two step. Thus, in the 
specific heat, we observe two peak temperatures. The higher 
and lower ones are defined as Tp2 and Tk, respectively. 

In Fig. 6(b), we plot Tp 2 and Tk as functions of A. Both 
temperatures are increased monotonically up to A = 0.5, but 
for A > 0.5, those seem to be almost constant, even if we 
change A. The higher peak temperature Tp 2 is found to be 
well fitted by AE3, which is the excitation energy in the three- 
site system, including one impurity and two conduction sites. 
The ground state of the three-site system is found to be char¬ 
acterized by ^tot = 7/2 for the case of A = 0. When we 
include the spin-orbit coupling, the octet ground state is split 
into three, two doublets and one quartet. Among them, one 
doublet becomes the ground state and the quartet is the first 
excited state. The energy difference between them is here de¬ 
fined as AE3. Again this is the level splitting due to the effect 
of hybridization. Note that the A dependence of AE2 is simi¬ 
lar to that of AE3, but the magnitude is found to be different 
fromTp2. 

As we observed in Fig. 6(b), the higher peak Tp2 is scaled 
by AE3. The lower peak is considered to be the Kondo tem¬ 
perature Tk. We find that Tk is increased monotonically up 
to A = 0.5 and it becomes constant for A > 0.5. In the limit 
of A = 00 , the j-j coupling scheme becomes exact and we 
know that one electron is accommodated in j = 7/2 octet. In 
such a limit, as emphasized in the section of introduction, the 
hybridization of one j = 7/2 electron with the conduction 
band leads to the Kondo effect. For enough large A, we easily 
deduce that the situation is not so changed from that in the 
limit of A = 00 and thus, it seems natural that the present Tk 
is constant in the region of large A. 

However, it is surprising that the same Tk as in the j-j cou¬ 
pling limit is obtained even for A = 0.5, which is not consid¬ 
ered to be large enough. When we decrease A, Tk is decreased 
rapidly, but we still find Tk in the present temperature range 
for A = 0.1. Thus, we arrive at an important conclusion that 
the picture of the Kondo effect on the basis of the j-j cou¬ 
pling scheme is applicable even for the realistic situation with 
the spin-orbit coupling in the order of X/U = 0.1 for elec¬ 
tron systems. The component of the j-j coupling scheme in 
the many-electron wave function is more persevering than we 
have naively expected. 

4.3 Effect of CEF potential 

Thus far, we have considered the situation without the CEF 
potential. As mentioned above, due to the effect of the hy¬ 
bridization, the level splitting has been found to occur. Next 
we include explicitly the CEF potentials. In Figs. 7, we sum¬ 
marize the results of entropy and specific heat for three CEF 
ground states, which are controlled by x. As easily understood 
from Fig. 3(b), the cases of a; = 1.0, 0.0, and —1.0 correspond 
to the ground states of F/", F^, and F^ , respectively. Since V 
is chosen to be very small here, we find the localized CEF 
states in the present temperature range. The residual entropy 



temperature T 


Fig. 7. (Color online) Entropy Simp (solid square) and specific heat Cimp 
(solid triangle) for (a) x = 1.0, (b) x = 0.0, and (c) x = —1.0. Other 
parameters are set as A = 0.2, U = 1, W = and V = 0.05. 


The entropy release from In 8 to In 4 or In 2 occurs and we 
find the peak in the specific heat at the corresponding temper¬ 
ature. The height of the peak depends on the released entropy, 
but the position of the peak is almost the same for three cases. 
The peak position is around 10“^*^ and this scale is determined 
by the CEF level splitting. As understood from Fig. 3(b), for 
A = 0.2, U = 1, and W = —0.001, the order of the CEF 
level splitting is in the order of 10“^. In the present calcu¬ 
lation, we set W = -10“"^, and thus, the order of the CEF 
level splitting is smaller in one order. In this sense, it is quite 
natural that the peak position appears around at T = 10“^ in 
common among Figs. 7(a)-7(c). 

Next we consider the multipole properties. After we faith¬ 
fully diagonalize the multipole susceptibility matrix eq. (17), 
we obtain the eigenvalues and eigenvectors. The multipole 
state with the largest eigenvalue is considered to be realized 
and the corresponding eigenvectors denote the optimized mul¬ 
tipole state given by the mixture of multipole components 
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Fig. 8. (Color online) Multipole susceptibilities for (a) a: = 1.0, (b) x = 
0.0, and (c) x = —1.0. Other parameters are set as A = 0.2, U = 1, 
W = -10-"^, and V = 0.05. 


with different ranks in the same symmetry group. 

First we consider the case without the CEF potential, i.e., 
W = 0. The results are not shown here, but the optimized 
multipole is always dipole for A = 0, since J is purely equal 
to total spin moment S = 7/2. When we increase the value 
of A, we observe the increase of other multipole components 
and in the j-j coupling limit, one electron in the /=7/2 octet 
carries all multipoles higher than dipole. Note that in a peri¬ 
odic system, the optimized multipole is affected by the lattice 
structure, electron hopping amplitude, and interactions. 

Next we include the CEF potential. In Figs. 8, we show the 
eigenvalues of multipole susceptibility matrix in each NRG 
step for the same parameters in Figs. 7. When we see the en¬ 
tropy and specific heat, it is difficult to find the difference 
in the cases for F/" and Fy ground states, corresponding to 
Figs. 7(a) and 7(c), respectively. In the plateau of In 8, there 
are no significant differences among three cases, but around 
at the temperature of the entropy release, it is possible to de¬ 


tect the different multipole states for three cases. For the case 
of F^ ground state, since F^ includes orbital degrees of free¬ 
dom, it carries higher-rank multipoles. On the other hand, Fg 
and F/ states carry only 4u multipoles, mainly characterized 
by dipole. 

However, for F/" and Fy ground states, we can observe the 
difference in the multipole state with the second largest eigen¬ 
value. Namely, they are 3g quadrupole and 5g quadrupole, 
respectively, for Figs. 8 (a) and 8 (c). The multipole with the 
largest eigenvalue is always 4u dipole, but the Curie constant 
for the multipole susceptibility is slightly different. Namely, 
in the j-j coupling limit, we obtain the dipole moment as 
(4/3)^b and (12/7)^b for Fg and F/ ground states, respec¬ 
tively. This difference appears in the values of Ex at the low¬ 
est temperatures in the present calculations. 

Thus far, we have considered the situation with small V 
such as E = 0.05, in order to focus on the localized proper¬ 
ties of /^-electron systems. Next we consider the case with 
lager V for the purpose to visualize the Kondo phenomenon. 
When we increase the values of V for Figs. 7(a)-7(c), we no¬ 
tice that the Kondo effect occurs for the case of Fy ground 
state, while for the cases of F/" and F^ ground states, we 
cannot observe the Kondo behavior even if we increase the 
magnitude of V up to 3.0. We remark that the Ou conduc¬ 
tion band is hybridized only with F/ state, since it includes 
(/i=2,a - /i=-2.a)|0)/v^> while Fg has no component of 
/m =2 cr|f*) fm=-2 cr)|f*) ^8 contains the Compo¬ 

nent of (/l=2,a+/l=-2,a) | 0 )/ Thus, in the following, 
we consider only the case of the Fy ground state. 

In Fig. 9(a), we depict the temperature dependence of en¬ 
tropy and specific heat for A = 0.2, U = 1, W = —10“^, 
X = —1.0, and V = 0.6. We find the plateau of In 8 around 
at T = 10 “^ and the entropy is changed to In 2 around at 
T = 10“^, forming the peak at Tp 2 in the specific heat. Then, 
the entropy In 2 is eventually released around at T = 10 “^, 
leading to the Kondo temperature Tk. The behavior of the 
entropy and specific heat are essentially the same as that in 
Fig. 6 (a) without the CEF potential. This is nor surprising, 
since the level splitting due to the hybridization plays the 
same role as that due to the CEF potentials. 

In Fig. 9(b), we plot the multipole susceptibility for the 
same parameters in Fig. 9(a). For T > Tk, we find qualita¬ 
tively the same behaviors as those in Fig. 8 (c). Around at the 
Kondo temperature, the Curie constant for 4u dipole gradually 
becomes zero, suggesting the screening of 4u dipole moment 
by the conduction electrons. Namely, the standard Kondo ef¬ 
fect is considered to occur in this case. 

When we change the values of V, we plot the peak temper¬ 
ature Tp 2 and the Kondo temperature Tk in Fig. 9(c). Again 
it is found that Tp 2 is well scaled by AT 3 , but we note that 
AT 3 converges to the local CEF excitation energy. Note here 
that the temperature is discrete in the NRG calculation, since 
it is defined as T = where A is the cut-off and 

N is the NRG step. Thus, the peak temperature Tp 2 includes 
an error-bar in the order of v/A. Due to the above reason, for 
small V, the solid circles seem to scatter around the horizontal 
broken line. 

However, for V > 1.0, Tp 2 is clearly larger than the CEF 
excitation energy and it seems to be scaled by AT 3 . For 
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Fig. 9. (Color online) (a) Entropy ^imp (solid square) and specific heat 
Cimp (solid triangle) for A = 0.2, U = 1, W = —10“^, x = —1.0, 
and V = 0.6. (b) Multipole susceptibilities for the same parameters in (a), 
(c) The peak temperature Tp 2 (solid circle) and the Kondo temperature Tk 
( solid squai'e) as functions of V. The solid curve denotes AE^, whereas the 
horizontal broken line indicates the local CEF excitation energy. The dotted 
cuiwe denotes the curve of IiiTk = b — alV"^ with appropriate constants a 
and b. 


and AEs is almost determined by the CEF excitation energy. 
When V is increased, the hybridization effect on the level 
splitting overcome the CEF potentials and AE 3 begins to de¬ 
viate from the local CEF excitation energy. We emphasize that 
Tp 2 correctly follow the above behavior of AE 3 . 

Finally, we discuss the Kondo temperature. As shown 
by the dotted curve in Fig. 9(c), Tk is well fitted by 
hexp(—a/V^) with appropriate constants a and b. This fact 
suggests that Tk actually indicates the Kondo temperature. 
Note that for V > 1.0, Tk becomes so high that we cannot 
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5. Discussion and Summary 

In this paper, we have discussed the Kondo effect in 
electron systems on the basis of the seven-orbital Anderson 
model by using the NRG technique. Note that we have con¬ 
sidered the single Ou conduction band. We have clarified that 
our understanding on the Kondo effect in the Ce compound 
also works in the /^-electron system with the spin-orbit cou¬ 
pling. If we are simply based on the j-j couping scheme, the 
result seems to be almost obvious. Namely, in the limit of 
large spin-orbit coupling, we accommodate one electron in 
Z=7/2 octet, while j=5/2 sextet is fully occupied. When one 
electron in the octet is hybridized with conduction band, it is 
possible to understand the Kondo effect in a similar way as 
that for the Ce compounds. 

From the qualitative viewpoint, the j-j coupling scenario 
may work for the understanding of the Kondo effect in Eu 
compounds. However, it is unclear whether the scenario is re¬ 
ally valid for actual compounds with the finite value of the 
spin-orbit coupling. In particular, in actual rare-earth mate¬ 
rials, the effect of Coulomb interaction is generally stronger 
that that of the spin-orbit coupling. Thus, naively speaking, it 
is difficult to accept even qualitatively the scenario based on 
the j-j coupling scheme in the rare-earth compounds. It has 
been the motivation of the present paper to clarify this point. 
Our results have clearly suggested that the Ce-compound-like 
Kondo phenomena should occur for /^-electron systems in 
the region of X/U in the order of 0.1, which is the value of 
actual materials. 

Since the purpose of this paper has been to clarify the ap¬ 
pearance of the Kondo phenomena in /^-electron systems 
from the conceptual viewpoint, we have not discussed the 
quantitative explanation of actual Eu compounds. In order 
to promote further the study of the Kondo effect in Eu com¬ 
pounds, it is inevitable to include the effect of valence fluctua¬ 
tions, which has been ignored in this paper. It has been pointed 
out that the valence fluctuations are involved in the heavy- 
electron formation in Eu compounds.Interesting properties 
induced by critical valence fluctuations have been actively in¬ 
vestigated by Watanabe and Miyake.®^' 

A way to investigate the effect of valence fluctuations on 
the Kondo phenomena in /^-electron systems is to include 
further the inter-site Coulomb interaction between impurity 
and conduction sites. Although for the standard impurity An¬ 
derson model without orbital degeneracy, it is easy to include 
such an inter-site repulsion, in the present seven-orbital An¬ 
derson model, there are several possibilities to consider such 
orbital-dependent inter-site repulsion. Thus, the first step to 
promote the research of the effect of valence fluctuations on 
the Kondo phenomena is to construct the valid model includ¬ 
ing inter-site repulsions. It is one of future problems. 

We have considered one conduction band, but it can be 
validated, for instance, in cage-structure compounds such as 
filled skutterudites, since the band-structure calculations have 
revealed that the main conduction band composed of pnicto- 
gen p electrons possess Ou symmetry.®®^ If Eu ion becomes di¬ 
valent in the cage of filled skutterudites, it is interesting to in¬ 
vestigate Eu-based filled slitterudites. Another cage-structure 
can be also the candidate of the research of the Kondo effect 
in /^-electron systems. If it is possible to synthesize Eu-based 
1-2-20 compounds, it may be interesting. As for the theoreti- 
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cal task, it is necessary to increase the number of conduction 
bands. For instance, we consider conduction bands with 
double degeneracy. Since they are hybridized with states, 
it may be possible to obtain the two-channel Kondo effect in 
the Eu compounds. It is another interesting future problem. 

In summary, we have analyzed the seven-orbital Anderson 
model by employing the NRG technique. When we have in¬ 
cluded the spin-orbit coupling in the situation, we confirm the 
Kondo behavior similar to that in the Ce compound with in¬ 
creasing the magnitude of the spin-orbit coupling, as expected 
from the j-j coupling scheme. An important point is that the 
Kondo effect similar to the case of n = 1 occurs for realistic 
values of the spin-orbit coupling even in the case of n = 7. 
Even when we include the cubic CEE potentials, we have also 
found the Kondo behavior similar to Ce compounds for some 
CEE parameter regions. 
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